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Abstract 



O: 

ly-^ . For the non-conservative Caldirola-Kanai system, describing a quantum damped har- 

monic oscillator, a couple of constant-of-motion operators generating the Heisenberg-Weyl 
algebra can be found. The inclusion of the standard time evolution generator (which is 
(-H ■ not a symmetry) as a symmetry in this algebra, in a unitary manner, requires a non-trivial 

^ f extension of this basic algebra and hence of the physical system itself. Surprisingly, this 

extension leads directly to the so-called Bateman dual system, which now includes a new 
particle acting as an energy reservoir. In addition, the Caldirola-Kanai dissipative system 

■ can be retrieved by imposing constraints. The algebra of symmetries of the dual system 
is presented, as well as a quantization that implies, in particular, a first-order Schrodinger 
equation. As opposed to other approaches, where it is claimed that the spectrum of the 
Bateman Hamiltonian is complex and discrete, we obtain that it is real and continuous, 

^ ■ with infinite degeneracy in all regimes. 

■ PACS: 03.65.-W, 02.20.-a, 2.30.Hq. 

O' 

O ■ 1 Introduction 

■ The interest in dissipative systems at the quantum level has remained constant since the 
early days of Quantum Mechanics. Within its framework, the attempts to describe damping, 
which intuitively could be understood as a mesoscopic property, have motivated a huge amount 

^ I of papers, both from a theoretical and a more applied points of view. 

■ Applications of quantum dissipation include quantum optics [1] and the study of decoherence 
phenomena [2]. It is customary to model dissipation by means of the theory of open systems 
or the thermal bath approach, in which a damped system is considered to be a subsystem 
of a bigger one with many, or infinite, degrees of freedom [3]. However, damped systems are 
interesting in themselves as fundamental ones if one wants to capture the essential dissipative 
behavior in a simple model. In particular, the quantum damped harmonic oscillator has attracted 
much attention, as could be considered one of the simplest and paradigmatic examples of non- 
conservative system. 

The quantum damped harmonic oscillator has been described frequently by the Caldirola- 
Kanai model [U [5], which includes a time-dependent Hamiltonian. It has been claimed that, 
in the quantum theory built out of this kind of models, uncertainty relations are not preserved 
under time evolution and could eventually be violated |6l [7] , but this apparent inconsistency is 
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associated with a confusion between canonical momentum and "physical" momentum (see, for 
instance, [8]; in fact, the same author proposed a logarithmic nonlinear Schrodinger equation to 
deal with this situation [9]). On the contrary, the main drawback of these models might be that 
they do not describe the interplay between the system and the environment, so that phenomena 
such as decoherence are not taken into account. It is nevertheless interesting to note that this 
kind of models can in fact be derived from a system-plus-reservoir model |10| [TT] in the case of 
Ohmic spectral density. 

The analysis of damping from the symmetry point of view has proved to be especially fruitful. 
In a purely classical context, the symmetries of the equation of the damped harmonic oscillator 
with time-dependent parameters were found in |12j . For the usual damped harmonic oscillator, 
finite-dimensional point symmetry groups for the corresponding Lagrangian (the un-extended 
Schrodinger group I14 | [TC]) and the equations of motion (S'L(3,M)) respectively, were found 
in |16| 117) . as well as an infinite contact one for the set of trajectories of the classical equation. 

As far as the quantum theory is concerned, in |18| the authors provided a neat framework 
to study a generalized version of the Caldirola-Kanai model with time-dependent parameters 
(frequency, damping coefficient and external force) , based on a quantum version of the Arnold 
transformation (QAT) |19) . The integrals of motion and symmetries were identified and exploited 
to calculate wave functions, basic operators and the exact time evolution operator, generalizing 
some of the results in |20) . Many papers related to the Caldirola-Kanai model keep appearing 
\21\ 122). showing that the debate about fundamental quantum damping is far from being closed 
(the reader can see |18| for further references on the Caldirola-Kanai model). For instance, at 
the quantum level of this model, it remains to address the symmetry role of the time translation 
generator ih-^. In fact, ih-^ acting on a solution of the Caldirola-Kanai Schrodinger equation 
is no longer a solution. As a consequence, the time evolution operator U does not constitute 
a one-parameter group of unitary transformations. Equivalently, the solution of the equation 
ih-^U = H{t)U is not e~T^^^''^\ but rather the time ordered product Te~i ^ ^^^^'^^ , referring to 

the Neumann series, or the Magnus series e~ft^''*^ |23) . 

The purpose of this article is to throw some light on the subject of quantum dissipation 
with the guide of symmetry. The starting point will be the Caldirola-Kanai system and we shall 
see that this model contains, somehow, a reminiscence of more general, system-plus-reservoir 
models. Using the QAT, there can be found basic integrals of the motion closing a Heisenberg- 
Weyl algebra, imported from those of the free particle, together with their squares. Time 
translations in the damped system do not belong to the imported quadratic, conserved operators. 
This is to be expected, as the classical equation of motion includes a friction term and the 
energy in this system is not conserved. The following question immediately arises: is there 
any finite-dimensional group of symmetries containing time translations and, at least, the basic 
operators? The answer is "yes", and we shall pay attention to this question in the case of the 
damped harmonic oscillator and the surprising consequences of the subsequent calculation: for 
this symmetry to act properly, it is necessary to enlarge the physical system with a new degree 
of freedom, corresponding to a new particle with interesting properties. 

In his original paper |24j , Bateman looked for a variational principle for equations of motion 
with a friction term linear in velocity, but he allowed the presence of extra equations. This trick 
effectively doubles the number of degrees of freedom, introducing a time-reversed version of the 
original damped harmonic oscillator, which acts as an energy reservoir and could be considered 
as an effective description of a thermal bath, as opposed to the standard approach followed in 
the theory of open quantum systems, where the thermal bath contains many, or infinite, degrees 
of freedom. The Hamiltonian that describes this system was later rediscovered by Feschbach and 
Tikochinsky |25| [26l [271 EH] and the corresponding quantum theory was immediately analyzed. 

Some issues regarding Bateman's dual system arose in the widespread Feschbach-Tikochinsky 
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construction of the corresponding Hilbert space and spectrum. This construction leads to a set 
of complex eigenvalues of the Hamiltonian (see |29) and references therein), and it has been 
interpreted that the vacuum of the theory decays with time |3oJ^ . In this paper we shall provide 
a rigorous construction of the physical Hilbert space and the real spectrum of the Hamiltonian, 
which are not found by means of the Feschbach-Tikochinsky construction. To do so, we shall take 
advantage of the adoption of a canonical quantization in which the Schrodinger equation is first- 
order. This approach is very similar to that presented in |31| and it also has some relationship 
with the one in |29) : there, the generalized eigenvectors corresponding to the complex eigenvalues 
of the Hamiltonian are interpreted as resonant states. 

Bateman's dual system is still discussed frequently |32| 133). Many authors have considered 
this model as a good starting point for the formulation of the quantum theory of dissipation. One 
of the aims of this paper will be to show that the study of the symmetries of the Caldirola-Kanai 
model leads to Bateman's dual system, so that it should be considered as a natural starting point 
for the study of quantum dissipation. This could be seen as a "step-back-and-forth" approach 
(from a model of bath, then the Caldirola-Kanai model and finally Bateman's dual system), 
arriving at an idealization which intends to capture some of the features of quantum dissipation. 
Moreover, we shall show that it is possible to get back to the Caldirola-Kanai system by imposing 
constraints in the Bateman classical system. 

The paper is organized as follows. We begin in Section [2] by recalling the results from |18| . 
in the case of the damped harmonic oscillator with constant coefficients in order to import basic 
operators from the free particle system, which satisfy the condition of being integrals of the 
motion and close a Heisenberg-Weyl algebra. In Section [3] we look for a minimal algebra of 
operators which includes the basic conserved operators of the Caldirola-Kanai system and the 
generator of time translations ^. In so doing, we arrive at an algebra of conserved operators of 
an enlarged physical system, which turns out to be the Bateman dual system. We also show how, 
at the classical level, the Caldirola-Kanai system can be recovered from the Bateman system by 
a constraint. 

In Section m we revise the quantization of the Bateman dual system. In particular, we find a 
first-order Schrodinger equation, from which the wave functions and the energy spectrum can be 
obtained. Finally, an Appendix is devoted to the study of a non-minimal, infinite-dimensional 
symmetry algebra for the damped particle. 



2 Basic operators in the Caldirola-Kanai model 

The classical equation of motion describing a Damped Harmonic Oscillator (DHO) is given 

by: 

X + ■yx + uj'^x = , (1) 

where a; is a frequency and 7 is the damping constant, both defining the system. This equation 
can be derived, in particular, from the Hamilton principle with a time-dependent Lagrangian. 
The corresponding, time-dependent, Hamiltonian function, to be named Hck, is given by 

HcK = e-^'^ + lmu^xh'yK (2) 
2m 2 

The canonical quantization of Hqk leads to the Caldirola-Kanai equation [HE], which is a 
Schrodinger equation for the DHO: 

^^9^2^^ "2^^ ^ 2"^""^ 

^We feel that the ultimate reason of this interpretation is the lack of a vacuum representation of the relevant 
symmetry algebra (see Section [3}. 
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Some comments on the features of this equation have aheady been made in the Introduction. 
Due to the explicit time dependence of the Hamiltonian function, the corresponding operator 
HcK is not a conserved operator and it is not possible to formulate an eigenvalue equation for it 
(that is, a time-independent Schrodinger equation). Hence, there does not exist a one-parameter 
group of time evolution. This is directly linked with the fact that the operator HcK does not 
preserve the space of solutions of the equation Q (see |18| for a further discussion) . 

Now we are just interested in identifying basic position and momentum operators associ- 
ated with classical conserved quantities (Noether invariants). In general, a conserved quantum 
operator 0{t) must satisfy the relation: 

^d{t)^^0{t)+'-[H{t),d{t)]=Q. (4) 

These operators are also known as dynamical invariants. Although the particular method to 
obtain these conserved operators is not fundamental in what follows, we find especially useful 
for that purpose the Quantum Arnold Transformation (QAT) technique developed in |18) . 

The QAT for the DHO maps unitarily solutions of the Schrodinger equation for the Galilean 
free particle to solutions of ([3]) and, accordingly, conserved operators in the Hilbert space of 
solutions of the free Schrodinger equation to operators acting on the Hilbert space of solutions 
of ©. 

In particular, we are going to map the Galilean momentum operator vr, corresponding to 
the classical conserved quantity 'momentum', and the position operator k, corresponding to 
the classical, conserved quantity 'initial position', to operators for the Caldirola-Kanai system. 
These free operators are, explicitly: 

d ^ ih d , , 

vr = -th— , K = K-\ r— , (5) 

OK m OK 

that is, those basic, canonically commuting operators with constant expectation values for the 
free evolution. 

The only relevant part of the QAT in which we are interested here is how the above operators 
are transformed (we refer the reader to [18] for further information). The QAT is defined just 
by two independent, classical solutions of ([1]), ui{t) and U2{t), satisfying the initial conditions 
ui(0) = 0,U2(0) = l,Mi(0) = 1,U2(0) = 0. Their Wronskian W{t) = ui{t)u2{t) - ui(t)u2{t) 
does not vanish. These classical solutions are: 

lii(t) = -^e'a* sinOt, ti2(t) = e"2*cosOt + -^e"2*sinQt, (6) 
for which W{t) = e"'''*, and 



is the characteristic frequency for the damped oscillator. Note that these solutions have good 
limit in the case of critical damping uj = (i.e. = 0). 

Performing the QAT for the operators ^ results in the formulas (see 



p = -ifiu^-^ - X , X = ^x + ih—-^ . (8) 
ox W W m Ox 

It can be checked by direct computation that these expressions satisfy the condition ^ provided 

that ui and U2 are solutions of ([TJ. Explicitly: 

P = -ih^—^ (2fi cos + 7 sin $7t) — + (7^ + 40^) sin Qt x , (9) 

2U ox 4il ' 

X = %-(2n cosOt-7 sinVLt) X + ih^—- smVLt— , (10) 
2il mil ox 
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with 



X,P =ih. (11) 



They are the basic, canonically commuting operators with constant expectation values (or, in 
other words, the two independent dynamical invariants) for the evolution of the Caldirola-Kanai 
system describing the quantum DHO. 

This way, we have imported the Heisenberg-Weyl algebra of conserved operators, which are 
symmetry generators, from the free particle system, but we have not done so for a generator of 
time evolution. The time evolution generator present in the free particle, and which is a constant- 
of-motion operator, may also be brought into a constant-of-motion operator in the Caldirola- 
Kanai system and will be realized as a quadratic operator in P and X (more precisely, P^), but 
this does not represent the actual time evolution generator of the Caldirola-Kanai oscillator. It 
is desirable that the system be completely characterized by its symmetries, including its time 
evolution, but so far, this had not been the case. 



3 Addressing dissipative systems in a symmetry framework 

Even though it is possible to set up a clear framework to deal with the quantum Caldirola- 
Kanai system for the DHO by employing the QAT, this does not provide by itself a well-defined 
operator (on solutions) generating the actual time evolution. As mentioned previously, this is 
rooted in the fact that the conventional time evolution generator is not included in the symmetry 
algebra: the Hamiltonian does not preserve the Hilbert space of solutions of the Caldirola-Kanai 
Schrodinger equation. 

We would like to emphasize that the underlying motivation for constituting a closed algebra of 
quantum observables is the general consensus that the quantization of any physical system must 
be a unitary and irreducible representation of some Poisson sub-algebra of classical observables 
characterizing the system to be quantized, irrespective of the particular method devoted to the 
achievement of this task (quantization method). The inclusion of the Hamiltonian in the chosen 
Poisson sub-algebra, and therefore in the quantum algebra, guarantees the characterization 
of the physical system. One may accordingly wonder what happens if the generator of time 
evolution for the damped harmonic oscillator is forced to belong to the algebra containing the 
basic operators. We shall pursue this issue in this Section. 

3.1 Including time symmetry 

In the Caldirola-Kanai model of the damped harmonic oscillator, neither the operator ih-^, 
nor HcK (which coincides with the former on solutions) preserve the space of solutions of the 
Schrodinger equation. We shall impose the condition of ih-^ being a symmetry generator. But 
it will be done in an elegant way, trying to close with X and P an enlarged Lie algebra of 
(conserved, symmetry generating) observables acting on the (possibly enlarged) Hilbert space 
Ti. Implicitly, we are forcing the system to become conservative, but this is the only information 
we are going to provide. 

As a first step, we shall take advantage of the explicit expressions obtained for the conserved 
basic operators: we compute the commutators of X and P with H = ih-^, and the resulting 
expressions are considered as new operators. After that, the commutators of the new operators 
with the formers are also computed, in the hope that this process ends up closing a Lie algebra 
at a finite number of steps. In fact, the resulting enlarged algebra is finite-dimensional and 
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includes X,P,H and four more generator^ (plus the central one /), denoted by Q,Il,Gi and 
G2, 

_7t , ^ d uj^ 21 

P = —the 2 (cos fit H — — sin Qt) — — h m — — e 2 sin fit x 
^ 2fl ' dx fl 

_ 2i 

- It 1 ■ \ -^^ ^ ■ d 

A = e 2 (cos — sm iZt) x + in — — sm iZt-— 

^ 2n ^ mfl dx 

_Tt_ , 7 , (9 bp' it_ 

H = —the 2 (cos fit — sin ritj— — \- m— e ^ sinQtx 

L CJ X ^ L 

~ rt, 87 6 2 

Q = e 2 (cos i7t — ■ sin fit) x + i?i — — sin fit—— 

2i2 mu ox 



Gi 
G2 



1 



~ 4172 + 7^ cos 2fLt + 27O sin 2J7t) 



02 '^^"^ 



which close the eight-dimensional algebra A: 



X,P 
H, Gi 



-Go 



ihl 

ih 

— ( 

m 

—ihGi + ih'yG2 



-n 



m 



ih 



-2ih'yX P + ih-yQ 

m 

-ih^Gi + 2inuPG2 



Q,n 
x,n 
An 

H,P 
H, G2 



2ihGi - ihl 
ihGi 

—ihmupG2 

2ihmuP'X — ihmuj'^Q (12) 
2)ihmuPX — 2ihmuj'^Q — ih'jll 
-2ihGi + ih'yG2 + 2ihi . 



The outcome of this process depends, obviously, on the way ^ acts on the explicit realization 
of the Heisenberg-Weyl generators, which in this case encodes somehow the time evolution of 
the Caldirola-Kanai oscillator {X and P were conserved operators under the evolution of this 
system). 

Now we have to unveil the physical content of (|12|) . This algebra A corresponds to a central 
extension of an algebra, to be named A, and this particular central extension itself determines 

the actual basic conjugated pairs, fixing a specific quantization. The operators Q and 11 (plus 
I) expand a Heisenberg-Weyl subalgebra, and H, Gi and G2 expand a 2-D affine algebra (with 

H acting as dilations). However, in this realization Q and H are not basic operator^, while 
H and G2 are conjugate ones: time proves to be a basic variable, this way contributing to the 
parametrization of the solution manifold. This result is puzzling: time is not expected to be a 
coordinate or momentum of any degree of freedom and therefore the corresponding generator 
should not appear as an element of a conjugate pair. We then wonder whether it is possible or 
not to take advantage of the information encoded in the algebra A to describe a physical system 
in which we have an evolution with respect to the ordinary time variable. It could not even be 

discarded the possibility that in the quantum reduction process, in which Q and 11 cease to be 



^At least in the simpler case of the damped particle, infinitely many new generators can be included in its Lie 
algebra. See Appendix |X] for further details. 
''This can be seen as an anomaly, see [341 135] . 
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basic operators, H would regain its expected role of time evolution generator. However, this 
requires a deeper analysis and fortunately this is not the only solution. 

In fact, A is not the only possible central extension of A and our strategy here will be to 
consider other possible quantizations of the un-extended algebra A, labeled by its possible central 
extensions. A detailed study shows that there are three relevant kinds of central extensions, 
describing systems with different degrees of freedom. Thinking of the algebra above as an 
abstract Lie algebra, it can be shown that a couple of parameters, in addition to the ordinary 
mass parameter m, control the central extensions that are allowed by the Jacobi identity of Lie 
algebras. We shall focus on a one-paremeter family Ak which contains the most relevant cases 
from the physical point of view: 



H,X 
H, Gi 



ihl 

ih 



Go 



ni 



-ihGi + ih-fG2 + ih{k - 1)1 



m 



ih _ 



Q,n 
x,u 
An 

H,P 
H,U 

H, G2 



-2ih-fX P + ih'yQ 

m 

-ih^Gi + 2inu?G2 
It is now convenient to perform the linear shift: 

Q = Q + {k-l)X, 
so that the actual degrees of freedom diagonalize: 



2ihGi - ihkl 
ihGi 

—ihmuj^G2 
2ihmu!^X — ihmuj'^Q 
'iihmuP'X — 2ihmuj'^Q — ihjfl 
-2ihGi + ih-^G2 + ih{l + k)i . 



X,P 
X,Q 
P,Q 
H,X 
H,Q 

H, Gi 



ihl 
ih 

— ( 

m 

—ihGi + ih'yG2 

_^h^ 
m 

-ih-yil + k)X - 

ih 



-G2 



ih 



-p 



m 



+ ihjQ + —(1 - k)U 
m 

-ih-fGi + 2ihio^G2 



x,n 
An 

H,P 
H,U 

H, G2 



ih{k + l)Gi - ihkl 
ihGi 

—ihmuP'G2 

ihmuj'^{l + k)X — ihmuj'^Q 
ihmu?{2k + l)X 

— 2ihmoj^Q — ih'jll 
-2ihGi + ih-fG2 + ih{l + k)i . 



By noticing the appearance of the central generator / on the right hand side of these com- 
mutation relations, we can see that the centrally extended algebras Ak can be classified as 
follows: 

• For k 7^ ±1, a generic family Ak describes systems with 3 degrees of freedom: (X,P), 
(Q,n) and (H,G2), thus time being a conjugate variable. 

• For k = 1, Ai = A (already considered) describes an anomalous system with 2 degrees of 
freedom: (X,P) and {H,G2), time being again a conjugate variable. 
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• For k = —1, A^i describes a system with just 2 degrees of freedom: {X,P) and {Q,Il). 

The third case is what we are looking for: ahhough it contains two degrees of freedom, time 
is not a conjugate variable to any other. We should stress at this point that A-i is in fact the 
unique central extension, inside the family Ak, in which H is directly an ordinary time evolution 
generatoi0. This is the reason why the whole process is unambiguously defined. 

A complete study of the entire family of central extensions of A and specially the original 
case k = 1 and the corresponding aforementioned anomalous reduction of the quantum repre- 
sentation, which might eventually end up with a non-conserved Hamiltonian operator associated 
with H, will be published elsewhere. 

Explicitly, for A-^i we have: 



H,Gi 



-ih'jGi + 2ihh?G2 



X,P 


= ihl 


Q,U 


= ihl 


X,Q 


= Ct2 

m 


X,U 


= ihGi 


P,Q 


= —ihGi + ih'yG2 


An 


= —ihmijj^G2 


H,X 
H,Q 


m 

= -(-p + 2ri) 

m 
+ ih'jQ 


H,P 
H,Il 


= —ihmuJ^Q 

= -ihmu?{X + 2Q) 
— ih'ylil 



G2 



-2ihGi + ih'jG2 



In this case we notice that the operators Gi and G2 commute with the basic couples (X , P) 
and {Q, H) and, therefore, will be represented as constants (times the identity operator) in any 
irreducible representation. Even more, their commutation rules with H determine that the 
value of these constants must be zero {[H,Gi] = [-ff , G2] = 0, as Gi and G2 are proportional 
to the identity /; then, solve for Gi and G2 using the expressions abov^). Technically, the 
operators Gi and G2 are gauge, in the sense that, in the resulting physical system, their symmetry 
transformation does not produce any change in the parameters associated with basic operators. 
In consequence, the effective dimension of the algebra A-i is 5 -|- 1. Let us denote by B the 
effective reduced Lie algebra, whose homonymous generators {X,P), ((5,11), H and I have the 
commutation relations: 



X,P 


= ihl 


Q,n 


= ihl 


X,Q 


= 


X,Il 


= 


'q,p 


= 


p,n 


= 


H,X 


m 


H,P 


= —ihmuj'^Q 


h,q' 


= -{-P + 2U) 

m 

+ ihjQ 


H,U 


= -ihmu?{X + 2Q) 
— ih'fll . 



(13) 



^This is still true for the complete family of central extensions of the algebra A. In addition, for k = —1 both 
degrees of freedom share the same central extension parameters. 

^To be precise, in the critical-damping case Q, - 
still can be represented trivially. 



;2_ ^ 



this only implies Gi 



G2, although they 
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This way, we end up with two pairs of independent, canonical operators along with a Hamil- 
tonian. Since we have selected the representation in which H is not a basic operator, it can be 
written in terms of the basic ones in an irreducible representation (note that symmetrization 
has been imposed to solve the ordering ambiguity, so as to achieve unitarity): 



H = —UP - ^{QU + UQ) 

m 2 m 



2A2 



muj XQ — mijj Q 



(14) 



The physical content of this Hamiltonian can be analyzed classically by solving the Hamilton 
equations of motion. After that, a canonical quantization would also be possible. However, it is 
interesting to establish the relationship of the physical system at which we have arrived with a 
known one: the new system, with two degrees of freedom, is actually the Bateman dual system. 

3.2 Bateman's dual system 

Let us perform the following linear transformation of operators: 



X = y+ —Py + l-x 
m7 2 

F = Px — m—y — m — x 
2 7 

Q = -y Py + -X 

m7 2 

- 1 ^ ^"^ - 
II = Px + m—y — m — x . 

2 7 



(15) 



the inverse of which is given by: 



x = X + Q 
1 



Px 

y 



iP + U) + m—{X + Q) 
2 7 

1 



— (p-n) 

m7 

Py = P-tl+ im7(l - Q) . 



(16) 



This transforms the Hamiltonian (1141) into the so-called Bateman's dual Hamiltonian 



H = + ]r{ypy - xpx) + mO^xy . 
m 2 



(17) 



H closes a 5+1 dimensional algebra with {x,Px) and (y,py): 



[X,Px 



ihl 



= 

[y,Px] = 
ih 



H,x 
H,y 



J-py + mlx) 

ih . 7~N 
— {-Px - m-y) 
m 2 



[y,Py 


] =ihl 


[X,py 


]=0 


\Px,Py 


]=0 


H,px 


= ih{- 


H,Py 





(18) 
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It can be checked that the corresponding classical Hamiltonian H describes a damped particle 
{x,Px) and its time reversal {y,Py) by computing the second order classical equations out of the 
Hamilton equations: 

X + ■yx + uj'^x = , y — ■jy + Lo'^y = . (19) 

The system is conservative, so that our objective of including the generator of time evolution 
among the symmetries has been accomplished, though at the cost of including a new degree of 
freedom. 



3.3 Recovering Caldirola-Kanai system from the conservative Bateman sys- 
tem 

Before addressing the quantum description of Bateman's dual system in the following section, 
let us insist on the classical theory in this subsection. It is interesting to realize how the 
Bateman dual system actually generalizes the Caldirola-Kanai damped harmonic oscillator in 
the sense that a Caldirola-Kanai particle and its mirror image can be recovered by a canonical 
transformation. 

In Bateman's original paper |24| . the Caldirola-Kanai Lagrangian was obtained by imposing 
the simple constraint y = e^^x (Dekker [28] also described an analogous process by using a 
complex canonical tranformation) . Here we provide a (real) canonical transformation which takes 
the Bateman Hamiltonian to the difference of a couple of dual Caldirola-Kanai Hamiltonians. 
That is: 

/ 1 / 7t , 7 

y = ^[y-(^^ + TT-p^^ 

Py = ^[Py-^e ^'p. + m-x 
.' = ^(x + e-^,-^e-^ (20) 

/ ^ f ^'^ , -yt 7-7* 

P^ = ^ ( j^Px + py- m-ei x 
t' = t. 

This transformation, although explicitly time dependent, is the lifting of a canonical transfor- 
mation among initial constants. After performing the transformation, the new Hamiltonian H' 
(including the time derivative of the corresponding generating function) becomes: 

2m 2 2m 2 ^ ' ^ ' 

which is nothing other than the Caldirola-Kanai Hamiltonian in terms of the variables {x',p'^) 
minus a dual Caldirola-Kanai Hamiltonian in terms of the variables {y' ,p'y). The relative minus 
sign in the new Hamiltonian H' is unavoidable at least with a real canonical transformation. 

A constraint to obtain the original Caldirola-Kanai Hamiltonian now reduces simply to the 
form: 

y' = 0, p'=0. (22) 



4 Quantization of Bateman's dual system 

In principle, the quantum theory of the Bateman system can be addressed performing the 
usual canonical quantization in the position representation. Then, the Schrodinger equation for 
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the Bateman Hamiltonian (|17p is given by: 



- -'H -^-) • <-) 

It has been argued that the quantum Bateman system possesses inconsistencies, such as 
complex eigenvalues and non-normalizable eigenstates \25\ \26\ [271 128| . However, it should be 
noted that this observation is meaningless if the proper Hilbert space of the system is not 
specified, together with the precise way in which it is constructed. In this respect, Chruscihski 
& Jurkowski |29) showed that H has in fact a real, continuous spectrum (we shall provide a proof 
of this in Subsection l4.2p . and that the complex eigenvalues are associated with resonances, which 
in last instance are the responsible of dissipation. 



4.1 First-order Schrodinger equation 

In order to construct the quantum representation of the Bateman algebra (jlSp . we shall 
choose a mixed representation of position- momentum on complex functions of x and py: 

d d 
X = x, px = -'^^-Q^ ' y = ' Py =Py (24) 

With this choice to represent the basic operators, the Hamiltonian becomes: 

H = ihilx -^)^ + ihi^Py + m^l^x)^ + ihl , (25) 
2 m ox 2 opy 2 



that is, a first-order operator which is Hermitian with this particular ordering. Operators (j24p 
and ([25]) (together with the identity) are Hermitian and represent the algebra ([18]) on complex, 
square integrable functions (j){x,Py) with integration measure dxdpy. As a consequence, the 
time-dependent Schrodinger equation becomes also first-order: 

4.2 Solving the equation: eigenfunctions and spectrum 

Before proceeding with the solution of the time-independent Schrodinger equation, let us 
firstly note that, being the Schrodinger equation first-order, it is possible to find the general 
form of the solution: 

(j)(x,py, t) = 62*/ ( e2*(xcos($7t) — ^ sin($7t)), e2*(p„ cos($7t) -|- milx sin(ilt)) ) , 

where / is an arbitrary complex function of two variables. This expression can be written in a 
more convenient form in order to compare with the solutions of the time-independent equation, 
by defining 

1 1 

Z+ = X + l—p^Py , Z-=X- i—^Py , (27) 

mil mil 
which are complex and z_ = z*^ when € M (the underdamping case). Now, the solution is 
written: 

(p{x,py,t) = 62*/ (^e(2+*f^)*z+, 6(?-*™)z_) , (28) 

or even: 



<^{x,py,t) = 6i*/ { 6^*^"^, e^'z+z. ) . (29) 
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Since this system is autonomous, it makes sense to formulate an eigenvalue problem for the 
Hamiltonian. It was already mentioned in the Introduction that several attempts of construct- 
ing the eigenfunctions of the Bateman dual system can be found in the literature. The most 
widespread is Feschbach-Tikochinsky's |25| 126): they in fact looked for creation and annihilation 
operators to find the Hamiltonian eigenfunctions. This approach is appealing, as it unveils an 
SU (1,1) structure for the system |30) . although the fact that the Hamiltonian is not positive- 
definite leads to undesirable consequences: an unphysical complex, discrete spectrum along with 
the corresponding (non-normalizable, not even to a delta function) wave functions. It should be 
mentioned that many authors have nevertheless considered that this construction leads to the 
correct physical description of the quantum Bateman system \28\ I30| . 

Two sets of authors proposed independently a different approach to find the real part of the 
spectrum (and therefore, the physically meaningful one) in the underdamping case. On the one 
hand, Rideau, Anderson and Hebda |3T] performed a canonical quantization in a mixed position- 
momentum representation (compare (f28|) with their resutlts). On the other hand, Chruscihski 
and Jurkowski |29| made use of angular variables to perform a first-order quantization of the 
Hamiltonian, arriving at the correct spectrum, although it is not clear from their construction 
to what extent they can give a consistent construction of the basic operators (compare ()29p 
with their resutlts). We shall follow the first approach, although we shall give the corresponding 
expressions for every regime: underdamping, overdamping and critical damping. 

The general solution of the time-independent Schrodinger equation can be found in terms of 
the variables (f27|l . The eigenfunction of H with eigenvalue E is given in terms of an arbitrary 
function g of a single argument: 

(|)--.(.,..(|)-"). 

This means that the spectrum is infinitely degenerate. To break the degeneration, we proceed 
as usual: we select another operator which commutes with the Hamiltonian and find their 
simultaneous eigenfunctions. In fact, it is possible to split the Hamiltonian into two operators 
of that kind: 



\m ox opy / 



(30) 



^ = ^^l(^i-+Py^ + ^]^ (31) 



2 \ dx "^'dpy 
with: 

H = Hn + D, [Hn,H] = 0, [D,H]=0. 
Then, the simultaneous eigenfunctions are, simply: 

which satisfy: 

Hn (j)n,X = nhVL (j)n,\ , D (pn,X = A/l7 (pn^X ^ H (j)n,X = KnVl + A7) (t)n,X , 

where the allowed values of n and A depend on the regime. 
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Underdamping. In this case, Q is real; therefore, n must be real since the spectrum of Hq 
must be real. Furthermore, is a pure phase, twice the one of z_. So, n is necessarily an 

integer. This may be seen as a consequence of the fact that Hq represents a rotation in the 
X, py plane for > 0. 

Regarding the value of A, it must be real, given that 7 is real. No further constraints are 
required for (j)n^\ to be well-defined. This is consistent with D representing a dilation, with 
spectrum the real line. So, the spectrum of the Hamiltonian is: 

E = nhn + A/i7 . (33) 

That is, we obtain a spectrum which is continuous and countably, infinitely degenerated. These 
results coincide with those in |29) . although they are obtained here in a neater way: in |29) the 
authors quantize angular variables and hence the basic operator "multiply by the angle" is not 
defined. We have avoided this problem: our basic operators are just (p4 



Overdamping. Now is pure imaginary, c.f. ([7|). In order to keep the spectrum real, n must 
be pure imaginary as well. Let us then define 

n = —in , Q = —iQ . 

The variables (j27p are now real: 

1 1 

z+ = x + —~-Py , Z_ = X ^py , (34) 

mil mil 

and it can be checked that solutions (1321) can be written in the form: 



(I)n,x{z+,z_) = e 'py'l^z+z-j . (35) 

Given that arccoth( ™^^ ) is not a periodic function, h can take any real value. This is consistent 

with the fact that the action of Hfi is analogous to that of a space-time boost in this regime. 

For the value of A, the same considerations as before lead to the conclusion that it must be 
real. Therefore, the spectrum of the Hamiltonian is now real with two real indices: 

E = hhQ + Xh'y . (36) 

In other words, the spectrum in this case is continuous and uncountably, infinitely degenerated. 



Critical damping. In order to find the eigenfunctions and the spectrum in the critical damp- 
ing case, it is necessary to solve the eigensystem for the operators pO|) and (|3ip when 0, is zero. 
This leads to eigenfunctions of the form: 

c^lA^,Py) = e'—[pl)---'\ (37) 

satisfying (^0 = ^n=o): 

Ho4>l^ = kh^<t^U. DcPl, = Xh^<Pl, H4>i^ = h^{k + X)ct>l^, 

where k and A must be real, so that the spectrum 

E = hjik + A) (38) 
is again continuous and uncountably, infinitely degenerated. 

It is now clear that this quantization, which renders the Hamiltonian first-order, allows a 
straightforward construction of the Hilbert space of the quantum Bateman dual system as well 
as the corresponding physical spectrum, which is not complex but real, as it must be for a 
Hermitian Hamiltonian. 
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5 Conclusions and Outlook 



In this paper we have found, starting with the non-conservative Caldirola-Kanai model for 
the quantum damped harmonic osciUator, a larger system containing two degrees of freedom 
(the original system and its time reversal one) by imposing that the time evolution generator 
ih^ closes a symmetry algebra with the basic operators X and P. The system obtained is 
the Bateman dual system which, in some sense, can be considered as the simplest model for a 
system-reservoir approach to the damped harmonic oscillator. 

It should be stressed that the process followed is unambiguous, in the sense that the symmetry 
algebra of the Bateman dual system is the only extension of the Caldirola-Kanai symmetry 
algebra where the time generator is not associated with a canonical variable. It would be worth 
investigating, however, the other possible extensions of the Caldirola-Kanai symmetry algebra 
where this is not primarily the case, for both the two and the three degrees of freedom cases 
(see Section \3.1\i . 

As commented in the Introduction, the quantization of Bateman's dual system has usually 
been considered full of inconsistencies since Feschbach and Tikochinsky |25| \26\ [271 EH] con- 
structed its Hilbert space made of non-normalizable energy eigenstates with complex discrete 
eigenvalues. Many other authors insisted on this construction and suggested alternatives |30|l32j. 
but in this paper we have shown that, using a first-order Schrodinger equation, the quantization 
is well defined and the spectrum is real and continuous, with infinite degeneracy |3H I29|. 

The only pathology of the Bateman dual system is that the spectrum of the Hamiltonian 
is not bounded from below, but this is not a problem as long as the whole Bateman system is 
isolated (i.e. there is no external interaction). This problem can be easily solved by restricting 
the initial conditions for the two particles in such a way that the kinetic energy be positive 
|24j . This effectively describes a system with just one degree of freedom, as shown in Section 
13.31 bv imposing the constraint ([22]) after performing the canonical transformation (j20p . which 
is equivalent to a relation among initial constants. It would be interesting to investigate the 
resultant one-degree-of-freedom systems after imposing other relations among initial constants 
that give rise to positive kinetic energy. 

The Bateman dual system solves many of the problems that the Caldirola-Kanai model 
for damped harmonic oscillator suffers from. For instance, while the time evolution under the 
Liouville equation for the density matrix of the Caldirola-Kanai system does not give rise to de- 
coherence (one of the main points raised against the Caldirola-Kanai model), the Bateman dual 
system suffers from decoherence and classical correlations after the mirror degree of freedom is 
traced out |36| . Also, selecting appropriate dense subspaces of the Bateman Hilbert space, an 
effective evolution describing damped oscillations with complex discrete frequencies is obtained 
|31| . and these complex discrete frequencies, appearing as the eigenvalues of the Bateman dual 
Hamiltonian in the Feschbach-Tikochinsky quantization, turn out to be resonant energies corre- 
sponding to resonant (Gamow) states for the Bateman dual system |29) . Even more, an splitting 
of the Hilbert space can be made in such a way that in each subspace the originally unitary time 
evolution falls down to a semigroup, and therefore irreversible to an evolution [29] . In addition, 
the Caldirola-Kanai system can be derived from Bateman's dual system (see Section \3.3^ . 

To end this section, we shall mention that the result that the algebra A ([12]) is finite- 
dimensional depends crucially on the fact that the frequency w and the damping coefficient 7 
are constant. For time dependent frequency u}{t) and damping coefficient 7(t), the procedure 
of iteratively commuting the resulting new operators with all other previously computed will 
not close to a finite algebra (except, perhaps, for very special cases). In this case the resulting 
extended algebra will probably contain an infinite number of degrees of freedom, a situation that 
presumably fits better with the concept of thermal bath. A deeper study of this general case 
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would be extremely interesting, at least in a perturbative way. 



A Appendix: Infinite-dimensional symmetry in the "free" damped 
particle 

In this appendix, we turn our attention to the "free" damped particle as the simplest case of 
physical system subjected to a dissipative force and perform a similar analysis to that carried 
out in Subsection 13.11 for the damped harmonic oscillator, forcing the introduction of the time 
generator in the symmetry algebra. 

The basic operators for the damped particle (obtained as the a; — )■ of the damped harmonic 
oscillator, see eq. ([QlllOp) 

P = -in^, X = x + —il-e-^')^. (39) 
ox mj ox 

Renaming P = Pq, we introduce operators P„ and (n being a non- negative integer) 

d - d 

ot ot 

ox 

X = x + —{l-e'')^. (40) 
m7 ox 

It is interesting that they close an infinite-dimensional Lie algebra with non-null commuta- 
tors: 

[Hg, Pn] = -ih-ynPn-i [Hdp, Pn] = -ih-ynPn 

[Hg, X] = -i-Po [Hdp, X] = -i-Pi 

m m 

[X, Pn] = hYn [Hdp, Yn] = -ihjnYn 

[HG,Yn\ = -ihjnYn-i [HDP,HG]=ih^HG , (41) 

containing a (centrally extended) Galilei algebra as a subalgebra, for which Yq = i is the central 
generator. Even, the generators in the left column of (j4ip . for n = 0,1, also close a finite- 
dimensional subalgebra in which Hg is a dynamical generator dual to Yi. 

It is interesting to look at the algebra (j4ip as some sort of gauge algebra. In fact, we may 
start from the specific realization Hg,X , Pq,Yq of the centrally extended Galilei algebra and 
try to turn it into a gauge (local) algebra by allowing the multiplication by certain functions of 
time (positive integer powers of e~'^^) of the space and time generators, and accordingly of the 
central one. This can be seen as the simplest symmetry of the free particle which results as a 
consequence of requiring the gauge invariance under the generator e~^^HG- 

It would be worth investigating the possibility of repeating this brief analysis for the case of 
the symmetry of the harmonic oscillator. 
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